Snow stratigraphic heterogeneity within ground-based passive microwave radiometer footprints: implications for emission modeling

Two-dimensional measurements of snowpack properties (stratigraphic layering, density, grain size and temperature) were used as inputs to the multi-layer Helsinki University of Technology (HUT) microwave emission model at a centimeter-scale horizontal resolution, across a 4.5 m transect of ground-based passive microwave radiometer footprints near Churchill, Manitoba, Canada. Snowpack stratigraphy was complex (between six and eight layers) with only three layers extending continuously throughout the length of the transect. Distributions of one-dimensional simulations, accurately representing complex stratigraphic layering, were evaluated using measured brightness temperatures. Large biases (36 to 68 K) between simulated and measured brightness temperatures were minimized (-0.5 to 0.6 K), within measurement accuracy, through application of grain scaling factors (2.6 to 5.3) at different combinations of frequencies, polarizations and model extinction coefficients. Grain scaling factors compensated for uncertainty relating optical SSA to HUT effective grain size inputs and quantified relative differences in scattering and absorption properties of various extinction coefficients. The HUT model required accurate representation of ice lenses, particularly at horizontal polarization, and large grain scaling factors highlighted the need to consider microstructure beyond the size of individual grains. As variability of extinction coefficients was strongly influenced by the proportion of large (hoar) grains in a vertical profile, it is important to consider simulations from distributions of one-dimensional profiles rather than single profiles, especially in sub-Arctic snowpacks where stratigraphic variability can be high. Model sensitivity experiments suggested the level of error in field measurements and the new methodological framework used to apply them in a snow emission model were satisfactory. Layer amalgamation showed a three-layer representation of snowpack stratigraphy reduced the bias of a one-layer representation by about 50%

Spatio-temporal influence of tundra snow properties on Ku-band (17.2 GHz) backscatter

During the 2010/11 boreal winter, a distributed set of backscatter measurements was collected using a ground-based Ku-band (17.2 GHz) scatterometer system at 26 open tundra sites. A standard snow-sampling procedure was completed after each scan to evaluate local variability in snow layering, depth, density and water equivalent (SWE) within the scatterometer field of view. The shallow depths and large basal depth hoar encountered presented an opportunity to evaluate backscatter under a set of previously untested conditions. Strong Ku-band response was found with increasing snow depth and snow water equivalent (SWE). In particular, co-polarized vertical backscatter increased by 0.82 dB for every 1 cm increase in SWE (R2 = 0.62). While the result indicated strong potential for Ku-band retrieval of shallow snow properties, it did not characterize the influence of sub-scan variability. An enhanced snow-sampling procedure was introduced to generate detailed characterizations of stratigraphy within the scatterometer field of view using near-infrared photography along the length of a 5m trench. Changes in snow properties along the trench were used to discuss variations in the collocated backscatter response. A pair of contrasting observation sites was used to highlight uncertainties in backscatter response related to short length scale spatial variability in the observed tundra environment

Evaluation of Operation IceBridge quick-look snow depth estimates on sea ice

We evaluate Operation IceBridge (OIB) ‘quick-look’ (QL) snow depth on sea ice retrievals using in situ measurements taken over immobile first-year ice (FYI) and multi-year ice (MYI) during March of 2014. Good agreement was found over undeformed FYI (-4.5 cm mean bias) with reduced agreement over deformed FYI (-6.6 cm mean bias). Over MYI, the mean bias was -5.7 cm but 54% of retrievals were discarded by the OIB retrieval process as compared to only 10% over FYI. Footprint scale analysis revealed a root mean square error (RMSE) of 6.2 cm over undeformed FYI with RMSE of 10.5 cm and 17.5 cm in the more complex deformed FYI and MYI environments. Correlation analysis was used to demonstrate contrasting retrieval uncertainty associated with spatial aggregation and ice surface roughness

Detached Eddy Simulation on the Turbulent Flow in a Stirred Tank

A detached eddy simulation (DES), a large-eddy simulation (LES), and a k-ε-based Reynolds averaged Navier-Stokes (RANS) calculation on the single phase turbulent flow in a fully baffled stirred tank, agitated by a Rushton turbine is presented. The DES used here is based on the Spalart-Allmaras turbulence model solved on a grid containing about a million control volumes. The standard k-ε and LES were considered here for comparison purposes. Predictions of the impeller-angle-resolved and time-averaged turbulent flow have been evaluated and compared with data from laser doppler anemometry measurements. The effects of the turbulence model on the predictions of the mean velocity components and the turbulent kinetic energy are most pronounced in the (highly anisotropic) trailing vortex core region, with specifically DES performing well. The LES—that was performed on the same grid as the DES—appears to lack resolution in the boundary layers on the surface of the impeller. The findings suggest that DES provides a more accurate prediction of the features of the turbulent flows in a stirred tank compared with RANS-based models and at the same time alleviates resolution requirements of LES close to walls

Historical Northern Hemisphere snow cover trends and projected changes in the CMIP6 multi-model ensemble

International audienceThis paper presents an analysis of observed and simulated historical snow cover extent and snow mass, along with future snow cover projections from models participating in the World Climate Research Programme Coupled Model Intercomparison Project Phase 6 (CMIP6). Where appropriate, the CMIP6 output is compared to CMIP5 results in order to assess progress (or absence thereof) between successive model generations. An ensemble of six observation-based products is used to produce a new time series of historical Northern Hemisphere snow extent anomalies and trends; a subset of four of these products is used for snow mass. Trends in snow extent over 1981–2018 are negative in all months and exceed −50×103 km2 yr−1 during November, December, March, and May. Snow mass trends are approximately −5 Gt yr−1 or more for all months from December to May. Overall, the CMIP6 multi-model ensemble better represents the snow extent climatology over the 1981–2014 period for all months, correcting a low bias in CMIP5. Simulated snow extent and snow mass trends over the 1981–2014 period are stronger in CMIP6 than in CMIP5, although large inter-model spread remains in the simulated trends for both variables. There is a single linear relationship between projected spring snow extent and global surface air temperature (GSAT) changes, which is valid across all CMIP6 Shared Socioeconomic Pathways. This finding suggests that Northern Hemisphere spring snow extent will decrease by about 8 % relative to the 1995–2014 level per degree Celsius of GSAT increase. The sensitivity of snow to temperature forcing largely explains the absence of any climate change pathway dependency, similar to other fast-response components of the cryosphere such as sea ice and near-surface permafrost extent

Effect of snow microstructure variability on Ku-band radar snow water equivalent retrievals

Spatial variability in snowpack properties negatively impacts our capacity to make direct measurements of snow water equivalent (SWE) using satellites. A comprehensive data set of snow microstructure (94 profiles at 36 sites) and snow layer thickness (9000 vertical profiles across 9 trenches) collected over two winters at Trail Valley Creek, NWT, Canada, were applied in synthetic radiative transfer experiments. This allowed robust assessment of the impact of estimation accuracy of unknown snow microstructural characteristics on the viability of SWE retrievals. Depth hoar layer thickness varied over the shortest horizontal distances, controlled by subnivean vegetation and topography, while variability of total snowpack thickness approximated that of wind slab layers. Mean horizontal correlation lengths of layer thickness were sub-metre for all layers. Depth hoar was consistently ~30% of total depth, and with increasing total depth the proportion of wind slab increased at the expense of the decreasing surface snow layer. Distinct differences were evident between distributions of layer properties; a single median value represented density and specific surface area (SSA) of each layer well. Spatial variability in microstructure of depth hoar layers dominated SWE retrieval errors. A depth hoar SSA estimate of around 7% under the median value was needed to accurately retrieve SWE. In shallow snowpacks <0.6m, depth hoar SSA estimates of ±5-10% around the optimal retrieval SSA allowed SWE retrievals within a tolerance of ±30 mm. Where snowpacks were deeper than ~30cm, accurate values of representative SSA for depth hoar became critical as retrieval errors were exceeded if the median depth hoar SSA was applied

Radio-frequency interference mitigating hyperspectral L-band radiometer

Radio-frequency interference (RFI) can significantly contaminate the measured radiometric signal of current spaceborne L-band passive microwave radiometers. These spaceborne radiometers operate within the protected passive remote sensing and radio-astronomy frequency allocation of 1400–1427 MHz but nonetheless are still subjected to frequent RFI intrusions. We present a unique surface-based and airborne hyperspectral 385 channel, dual polarization, L-band Fourier transform, RFI-detecting radiometer designed with a frequency range from 1400 through  ≈  1550 MHz. The extended frequency range was intended to increase the likelihood of detecting adjacent RFI-free channels to increase the signal, and therefore the thermal resolution, of the radiometer instrument. The external instrument calibration uses three targets (sky, ambient, and warm), and validation from independent stability measurements shows a mean absolute error (MAE) of 1.0 K for ambient and warm targets and 1.5 K for sky. A simple but effective RFI removal method which exploits the large number of frequency channels is also described. This method separates the desired thermal emission from RFI intrusions and was evaluated with synthetic microwave spectra generated using a Monte Carlo approach and validated with surface-based and airborne experimental measurements

Impact of measured and simulated tundra snowpack properties on heat transfer

Snowpack microstructure controls the transfer of heat to, as well as the temperature of, the underlying soils. In situ measurements of snow and soil properties from four field campaigns during two winters (March and November 2018, January and March 2019) were compared to an ensemble of CLM5.0 (Community Land Model) simulations, at Trail Valley Creek, Northwest Territories, Canada. Snow micropenetrometer profiles allowed for snowpack density and thermal conductivity to be derived at higher vertical resolution (1.25 mm) and a larger sample size (n=1050) compared to traditional snowpit observations (3 cm vertical resolution; n=115). Comparing measurements with simulations shows CLM overestimated snow thermal conductivity by a factor of 3, leading to a cold bias in wintertime soil temperatures (RMSE=5.8 ∘C). Two different approaches were taken to reduce this bias: alternative parameterisations of snow thermal conductivity and the application of a correction factor. All the evaluated parameterisations of snow thermal conductivity improved simulations of wintertime soil temperatures, with that of Sturm et al. (1997) having the greatest impact (RMSE=2.5 ∘C). The required correction factor is strongly related to snow depth () and thus differs between the two snow seasons, limiting the applicability of such an approach. Improving simulated snow properties and the corresponding heat flux is important, as wintertime soil temperatures are an important control on subnivean soil respiration and hence impact Arctic winter carbon fluxes and budgets

Are vegetation influences on Arctic–boreal snow melt rates detectable across the Northern Hemisphere?

The timing and rate of northern high latitude spring snowmelt plays a critical role in surface albedo, hydrology, and soil carbon cycling. Ongoing changes in the abundance and distribution of trees and shrubs in tundra and boreal ecosystems can alter snowmelt via canopy impacts on surface energy partitioning. It is unclear whether vegetation-related processes observed at the ecosystem scale influence snowmelt patterns at regional or continental scales. We examined the influence of vegetation cover on snowmelt across the boreal and Arctic region across a ten-year reference period (2000–2009) using a blended snow water equivalent (SWE) data product and gridded estimates of surface temperature, tree cover, and land cover characterized by the dominant plant functional type. Snow melt rates were highest in locations with a late onset of melt, higher temperatures during the melt period, and higher maximum SWE before the onset of melt. After controlling for temperature, melt onset, and the maximum SWE, we found snow melt rates were highest in evergreen needleleaf forest, mixed boreal forest, and herbaceous tundra compared to deciduous needleleaf forest and deciduous shrub tundra. Tree canopy cover had little effect on snowmelt rate within each land cover type. While accounting for the influence of vegetative land cover type is necessary for predictive understanding of snowmelt rate variability across the Arctic–Boreal region. The relationships differed from observations at the ecosystem and catchment scales in other studies. Thus highlighting the importance of spatial scale in identifying snow-vegetation relationships